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Abstract 

We study the quark-hadron phase transition at finite temperature with zero net 
baryon density by the Nambu-Jona-Lasinio model for interacting quarks in uniform 
background temporal color gauge fields. At low temperatures, thermal quark excita- 
tions are suppressed by a small expectation value of the Polyakov loop which prohibits 
excitation modes besides color singlet states, on the other hand at high temperatures 
quarks and anti-quarks are allowed to be excited almost freely with the Polyakov loop 
close to unity. Mesonic excitations are incorporated into the correlation energy as ther- 
mal fiuctuations of auxiliary bosonic fields. We show that at low temperatures these 
mesonic excitations dominate the pressure, while at high temperatures quarks and glu- 
ons instead dominate the pressure, mesonic excitations being absorbed in continuum of 
quark anti-quark excitations. We also show that there are two types of mesonic excita- 
tions, collective modes and non- collective modes, and only collective modes contribute 
to pressure at low temperatures. 
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1. Introduction 



Studying the phase diagram of quantum chromo-dynamics {QCDj^'^^ is one of the most 
fundamental problems in modern nuclear physics. One expects that at low temperatures 
and low densities, the chiral symmetry is broken spontaneously, and matter consists of 
interacting hadrons, where colored quarks and gluons are confined in individual hadron, 
while at sufficiently high temperatures and/or high baryon densities, the chiral symmetry is 
restored and matter is expected to turn into a plasma of deconfined quarks, antiquarks and 
gluons. Such extreme states of matter, commonly called the quark-gluon plasma, have been 
intensively studied experimentally by ultrarelativistic collisions of heavy nuclei with RHIC 
at Brookhaven and now with LHC at CERN.Sl' 

All QCD phase transitions are expected to occur in a non-perturbative regime with a 
strong QCD coupling. Although lattice QCD simulations by the Monte Carlo integration 
method have been shown very powerful in studying non-perturbative regime of QCD, it is 
plagued with the notorious sign problem in applying the same method to describe the region 
at finite baryon densities. For the exploration of the entire structure of the QCD phase 
diagram, we still need to resort to a more phenomenological approach adopting an effective 
theory which embodies some essential features of QCD. 

The Nambu-Jona-Lasinio (NJL) model was originally designed to describe profound con- 
sequences of presumed underlying chiral symmetry in hadron dynamics in analogy to the 
BCS theory of superconductivity.'^ It has later been transcribed to describe effective quark 
dynamics incorporating the symmetry.!^ The model, however, has no mechanism of quark 
confinement: therefore, when applied in the mean field approximation to a finite tempera- 
ture system, it generates thermal excitations of quark quasiparticles even at low temperature, 
where quarks are to be confined in individual mesons and baryons. 

To remedy this problem, a phenomenological model was proposed by FukushimaP by 
making use of the Polyakov loop which is known to characterize confinement- deconfinement 
transition in strong coupling pure gauge theories.® The Polyakov loop was inserted as a 
sort of fugacity factor in the quark distribution function and suppresses the quark single 
particle distribution by phase cancellations at low temperatures, leaving behind excitations of 
triad of quark quasiparticles in color singlet configuration. The quark-liberating deconfining 
transition was described as a change of the expectation value of the Polyakov loop treated 
as an order parameter with effective Landau- Ginzburg type free energy functional. 

Fukushima's original model has been reformulated as a mean field theory in uniform 
background temporal gauge field and has been studied extensively by others.'^''^''^ This 
model, referred to as the PNJL model, allows to incorporate mesonic correlation energies 
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beyond mean field approximation.'^''^ Such mesonic correlation energy had been studied 
previously in the original NJL model,'^''^ in analogy to the classic example of the RPA 
correlation energy in degenerate electron plasma.'^ The model has been used to compute 
other properties of hadronic matter, such as susceptibility at finite chemical potential,'^ 
thermodynamics in the regime of complex chemical potentiaP^ and the effect of confinement 
in the study of QCD phase diagram.!^ 

In this paper we reexamine the model in detail at zero chemical potential. We show 
that equation of state of a meson gas can be derived explicitly at low temperatures by 
the method of auxiliary fields which physically express effective meson fields build-up as a 
quark- ant iquark bound states as in the original Nambu-Jona-Lasinio model. Similar results 
have been obtained earlier by the other groups by numerical computations of the mesonic 
correlation energies with slightly different mo dels. '^''^''^''^^ We will show explicitly that in 
the chiral limit the equation of state becomes that of a gas of mesons which are described 
as collective modes of excitations and evaluate the correction due to the composite nature 
of mesons as the individual excitations of the quark triad. It is found that the former domi- 
nates the equation of state at low temperature and the latter do not alter the free meson gas 
result significantly. As the temperature increases the mean field contribution from under- 
lying quark quasiparticle excitations become important due to the non-vanishing Polyakov 
loop expectation value and will eventually dominate the free energy in deconfining phase. 
Collective mesonic excitations melt into the continuum of quark anti-quark excitations and 
the non-collective mesonic correlation remains negligible compared to the mean-field result. 

The rest of this paper is organized as follows. In the next section, we present the NJL 
model in uniform color gauge fields. We introduce bosonic fields as auxiliary fields in order 
to describe meson condensates and mesonic correlations and change four point NJL interac- 
tions to Yukawa type interactions. In section 3, the mean field approximation is reviewed. 
Under the mean field approximation, the fluctuation of meson fields are neglected. Therefore 
pressure is written in terms of constant auxiliary meson fields and thermalized quark quasi- 
particles. However these quark quasiparticles are not actually excited at low temperatures 
because the small expectation value of the Polyakov loop prohibits excitations of colored 
states. We discuss this mechanism in section 4 by focussing on the modification of the quark 
distribution function. In this section, we also note a change of the effective degrees of free- 
dom between the low temperature phase and the high temperature phase. In section 5, we 
calculate the contribution of mesonic correlation to the equation of state and show that they 
are decomposed into collective meson excitations and non-collective individual excitations of 
quarks or quark-triads. Especially in the chiral limit, where the bare quark mass vanishes, 
the contribution of collective modes is separated from that of non-collective modes and gives 
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the pressure of gas of free pions and free sigma mesons. We present our numerical results 
in section 6. It is shown that the pressure of free pion gas is converted continuously to 
that of free quark-gluon gas as the temperature increases, exhibiting a cross-over transition 
from hadron gas to a quark-gluon plasma. We also examine in detail the role of collective 
and non-collective degrees of freedom in low temperature confining phase and how they 
change as temperature increases, both in the chiral limit with massless pions and in more 
realistic situation with non-zero pion mass, by studying the modification of the dispersion 
relations. It will be shown that the collective meson poles are absorbed into the continuum 
of quark-antiquark excitations at high temperatures and individual excitations do not con- 
tribute significantly to the pressure of the system at all temperatures. We summarize this 
work in section 7. Some details of our computations are given in appendices. A condensed 
version of this work has been reported in.'^ 

§2. NJL model in uniform color gauge fields 

We begin with an Euclidean path-integral expression of the partition function for the 
NJL model of quarks in external temporal color gauge field: 

dr / d-^xCt^jL{q,q,A4 



(2-1) 



ZiT,A,) = j[dq][dq]exp 

where the effective Lagrangian is given by 

£NJL(g, q, A) = qi^YD^ - mo)q + G [{qqf + {iq-f^rqf] (2-2) 

for two flavor light quarks, q = {u, d) in external (classical) temporal color gauge field, 
= d^+gAoSfj^^Q and rriQ is the bare quark mass which breaks the chiral symmetry explicitly. 
Here we have used a standard Minkovsky metric notation with the real time being replaced 
by the imaginary time r = ixo and the time component of the gauge field replaced by 
= zAq^ We note here that quark fields are the only dynamical variables which describes 
thermal excitations of the system, while the temporal component of SU(3) gauge fields 
A4 = ^\°'A1, A'' being the 3x3 Gell-Mann matrices, only plays a side role imposing 
constraints on color configurations of thermal quark excitations. In the following calculation, 
we take Ansatz of diagonal representations for A4 along with the earlier works.'^''^''^''^ 

This partition function, Eq.f l2-ll) . contains the forth power of quark fields so that we 
cannot perform the integration over the fermion fields. For solving this problem, we introduce 
four auxiliary bosonic fields (pi = (a, tt) coupled to quark densities (gg, iq'~i5Tq) by multiplying 
Z(T, A4) by a constant Gaussian integral: 



j [d0]exp Jo j ^ '^'^^^ ~ "^"''^ + (tt + 2Giri^Tqf) 



(2-3) 
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where periodic boundary conditions are applied for each of the auxihary fields in the imag- 
inary temporal direction: 0j(r) = 0j(O), and the integration measure is chosen so that 
the integral gives unity. This procedure, known as the Hubbard- Stratonovich transforma- 
tion,'23'Eil' converts the four point NJL quark interaction to Yukawa coupling of the bosonic 
fields (pi to corresponding four quark densities. 



Z(T,A4) = / [dq][dq][dct>]exv 

where 



/3 

dr I d-^xCcs{q,qA,M) 



(2-4) 



£efr(g, g, 0, Ai) = q [iYD^, + a + i'y^T ■ tt] g - ^((a - niof + vrf ). (2-5) 

Performing the Grassmann integral over the quark Dirac fields, the partition function be- 
comes 

Z{T,A^) = j[d<p\e-'^'^^^^^ (2-6) 

where the exponent is given by 

A,) = ^jyrj d'x {{a - mof + tt^) - Trln [f3 (i7^D^ + a + ij,T ■ tt)] (2-7) 

where the trace in the second term is taken over the arguments of quark fields, including 
the space-time coordinates with anti-periodic boundary condition in the imaginary time 
axis, Dirac gamma matrices and the isospin and color indices. For constant auxiliary meson 
fields (((T, tt) = constants), the trace can be expressed as a sum over the fermionic (quark) 
Matsubara frequencies (e^ = (2n + 1)itT) and integrals over the spatial momenta: 



d^p 



(2-8) 



leaving behind only the trace over the color indecies. Here we have taken flavor sum with 
Nf = 2 assuming degeneracy of up and down quarks (m„ = rrid = ttiq). 

We evaluate the functional integral over the mesonic auxiliary fields by the method of 
steepest descent. Let 0o = {ctq, ttq) give local minimum value of the integrand so that 

51 , , 

— =0. (2-9) 

Hereafter we choose ttq = 0, along with ctq = —Mo, which evidently satisfy stationary 
conditions for effective pion fields. Shifting the integration variables by, ^pi = (pi — (po, we 
expand the integral in a power series of </?: 

, 1 5^1 



" ' 2 5(pi5(pj 



(PiiPj + --- (2-10) 
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The first order term of </? becomes zero because of equation deciding local minimum point. 
Keeping only up to quadratic term in the expansion, 



Z(T, A4) ~ e"^o / [dip]exp 



1 S^I 



2 6(()i6(j)j 

The gaussian integrals can be easily performed, yielding 



<t>=4>o 



(2-11) 



Z{T, A4) = e-^" [Det -^1 = e~'°~'^^''^''^ (2-12) 

d(pid(pj' 

which gives the thermodynamic potential: 

«(T.^,)=t(/o + 1tVm1u^) (243) 

up to one-loop fluctuations of the collective meson fields, where the trace is to be performed 
over the space-time coordinates of the auxiliary meson fields which obey a periodic boundary 
condition in the imaginary time direction. 

The first term of Eq. fl2-13p represents the thermodynamic potential under the mean field 
approximation and the second term represents the contribution of mesonic correlation to the 
thermodynamic potential. These mesons are free because we kept the expansion of A4) 
up to quadratic term in Eq.f l2-10|) . If we wish to take into account mesonic interactions, we 
need to keep higher order terms of expansion of A4). In this case we cannot integrate the 
partition function because of terms more than third power of (f so that we have to calculate 
the partition function perturbatively. In this paper, we ignore interactions between mesons. 

§3. Mean-field approximation 

The thermodynamic potential in the mean field approximation, f2yiF{T, A4), or the cor- 
responding pressure pmf(T, A4), may be retrieved from the leading term Iq by the relation: 

ilupiT, A4) = TIo = -PmfV. 

The explicit form of the leading term Iq is given by 

lo = ^y^iMo - mof - 2VNj ^ J ^trjn [p' ((e„ - gA,f + p^ + M^)] . (3-1) 
This leads to 

Pmf(T, a,) = -^(Mo - mo)' + 2TNj ^ J ^trjn [(5' ((6„ - gA,f + + Mo')] 

(3-2) 
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where e„ is fermionic Matsuibara frequencies, en = (2n + l)7rT, the trace is to be performed 
over the 3x3 color matrix ^4. The sum is taken over the Matsubara frequencies. Evalu- 
ating the discrete sum over the Matsubara frequencies by the standard method of contour 
integrationj^S we find 

Pmf(T, a,) = Ap(Mo) + 2NjT I ^tr, [in (l + e-^(^-+^^^^)) + In (l + e'/^^^^-^^^))] 

(3-3) 

where 

/MF(iVfo) = Sx2NfJ j^^E, - —(Mo - moY (34) 
is the pressure exerted by the zero point motion of the quark quasiparticles with energy 



Ep = yp^ + Mq in the "Dirac sea". The factor 3 in the first term of (13 4p accounts for color 
and the factor 2 for spin degeneracies. The second term in (13 4p is the pressure due to the 
chiral condensate which shifts the quark mass from rriQ to Mq and hence reduces the Dirac 
sea pressure. This term may be understood as the subtraction of the double counting of the 
mean-field effect as in the Hartree approximation. Note that this "vacuum pressure" does 
not depend on the external background gauge potential due to the cancellation of the effects 
of the potential on particle and antiparticle. It does not depend on temperature directly, 
but indirectly through the temperature dependence of the effective quark mass Mq. The last 
term in (13-31) is equal to the entropy of the system carried by the thermally excitations of 
quasiparticles multiplied by the temperature. It may be further reduced to more familiar 
form by partial integration, 

Pmf(T, a,) = p'mAMo) + 2Nj J -0^^tT, [f{Ep + tgA,) + f{Ep - tgA,)] (3-5) 

where 

nE,±igA,) = ^^^^^^l^^^^^ (3-6) 

is the quark (anti-quark) quasiparticle distribution function in the external gauge field po- 
tential. Note the change of sign in the gauge potential for quark {igA4) and for anti-quark 
{—igAi). The second term in (l3-5p corresponds to the pressure exerted by quark and anti- 
quark quasiparticles in the external gauge potential; the factor 2 accounts for spin degrees 
of freedom. The first term in (l3-3p exists even in vacuum at zero temperature so that it will 
be removed when we readjust the zero of the pressure. The second term present the pressure 
exerted by thermally excited quark quasi-particles which carries entropy. This term should 
be also removed in the low temperature phase where quarks are confined in hadrons. This 
can be done by choosing a specific form of as will be described in the next section. 
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In the mean field approximation, the value of the quark quasiparticle mass Mq is de- 
termined by the condition that the term hnear in the fluctuation 50 vanishes owing to the 
stationary condition fl2-9l) which may be written explicitly, for ctq = Mq, 



Mo - mo = 8GNf ^ J 



d^p Mq 



i27Tr{e^-gA,y+p^ + M, 



2 



(3-7) 



which is reduced to 

(1 - f{Ep + tgA,) - f{Ep - zgA,)) . (3-8) 
§4. Quenched quark distributions 

In the above calculations, the constant temporal gauge field A^ appears as a phase factor 
together with the quark quasiparticle energy in the single particle distribution function fl3-6l) . 
It looks very similar to the gauge invariant Polyakov loop phase integral. 



L(r) = Pexp 



■/3 

ig I drAi{r,T) 



(4-1) 



whose thermal expectation value measures an extra free energy associated with the color 
charge in fundamental representation fixed at a spatial point Although this connec- 

tion was first emphasized in the construction of the model by Fukushima, it looses its strict 
meaning in the present method, which follows the work of TUM group, since we assume that 
quarks are moving in a uniform background gauge field not fluctuating either in space or in 
imaginary time. Nevertheless we adopt the procedure of Fukushima, along with the TUM 
group, to replace the phase factor in the quark quasiparticle distribution function by the ther- 
mal average of the Polyakov loop. Taking diagonal representation L = (^g'/'i^ 6*"^^, e"**-'^^"'"'^^''), 
one find^ 

{-trJiE, + ^gA,)) -> ^ J f ^ ^ U{E,) (4-2) 

where 

^ = ^(tr,L) = ^((e*^^ + e^^^ + e-^(^^+^^))), i> = ^(tr.Lt) . (4-3) 

The anti-quark distribution, f^{Ep) = (|trc/(-Ep — 2(7^4)), is obtained from the quark distri- 
bution by the replacement of ^ by ^ and vice versa. At zero chemical potential, two distri- 
butions are identical and $ is real so that (l> = 3. We note this replacement corresponds to a 
Gaussian approximation in the statistical average, ignoring statistical correlations between 
the Poliakov loops, in the same spirit in the mean-field approximation. 
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Very interesting observation here is that although in the deconfining phase where ^ = 
(p = 1, one would simply have ordinary quark distribution, 

ME,)\,^, = ^^^^ (4-4) 
in the confining phase, where ^ = ^ = 0, we have instead, 

ME,)\,^o = ^^^^ (4-5) 

which may be interpreted as triad of three quark quasiparticles exciting together in color 
singlet configuration. It looks something like a baryon although no effect of interaction is 
taken into account between three quarks to form a baryon, as reflected by the same value 
of momenta for all three quarks. We note that the degrees of freedom of excitations is not 
reduced to 1/3 of color triplet free quark excitations. However, there still exist reduction 
of effective degrees of freedom due to the change of the phase space of quark momentum 
p and that of quark triad momentum 3p: The "triad" carries the energy E'tri = 3-Ep, the 
momentum ptri = 3p and the mass Mtri = 3Mq. This implies that the Lorentz invariant 
phase space integral over the quark momentum is reduced when written in term of the triad 
momentum, e. g. 

d^p 1 f d^Ptri 



Ep 9 J Etri ^ ^ 

Here we determine ^ phenomenologically by adding an effective potential of 

Q{T, <P) = {n{T, A^)) + U{T, <P) (4-7) 

where we deflne, following,'^ 

W(T, <P)/T^ = -]-h2{T)§<P - ^63(^>3 + + -M^^f (4-8) 
2 4 

with 

h2 (T) = ao + a,[^)+a2 (^) ' + (^) ' (4-9) 

The parameters of the effective potential are determined so that it gives prescribed values 
for ^ at given temperature in analogy to the free energy of the order parameter in the 
Ginzburg-Landau theory of superconductivity. The temperature dependent parameters in 
the above expression of U[$] are chosen so that <P = for the low temperature conflning 
regime where thermal quark quasi-particle excitations are forbidden, while ^ ^ 1 above 
deconflning temperature, quite opposite to the normal behavior of the order parameter: 
therefore the physical meaning of <P becomes clearer if it is called a " disorder parameter" . 
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We note here that, although gluon fields are not treated as dynamical variables, gluon 
excitations may be included in a phenomenological fashion through the effective potential 
by setting, 

W(r, ^ = 1) = (-ao - + h,] = -^T^ (4-10) 



3 " 4 7 90 

so that the equation of state matches to that of a free quark-gluon plasma at asymptotic 
high temperatures. 

With all these prescriptions, the equation of state in the mean-field approximation of the 
PNJL model is given as a maximum of 

PMF(T,^,Mo)=AF(Mo)-Z\pvac + 4x3x J -0^^f^{E,)-l{{T,<P) (4-11) 

with respect to variation of ^ and the quark quasiparticle mass Mq, which measures the 
magnitude of chiral condensate, the latter procedure gives the gap equation 



Mo - mo = 12GN 



which determines the optimal value of Mq. The constant Ap^^c is chosen so that the pressure 
becomes zero at zero temperature. In the confining phase, this equation could be interpreted 
as a gap equation for "baryon" (quark triad) mass, 

M,„ - 3mo = 3 X 12GiV^ x ^ J - f{E,„) - f{E,„)). (4-13) 

If we use a rescaled coupling G' = G/3, this equation indeed looks more like original Nambu- 
Jona-Lasinio gap equation for baryons. This interpretation however meets difficulty when we 
try to interpret the quark quasiparticle pressure fl4-lip as a pressure from thermal excitations 
of baryon and anti-baryons. By rescaling the momentum of quarks, the second term of fl4-lip 
may be rewritten in the limit <P = as 

4 X 3 X ^ [ ^^ Al. f^Et„) (4-14) 

The reduction of the degrees of freedom is by the factor 1/27, more than what we need for 
the compensation of the remaining color factor 3. 

We show the temperature dependence of ^ and Mq obtained from these mean-field equa- 
tions in Fig. 1. In the chiral limit mo = 0, the "gap equation" fl4-12p possesses a non-trivial 
solution Mq 7^ only at temperatures below a critical temperature Tc. In the ordinary NJL 
model, with no "Polyakov loop prescription", the chiral transition takes place at relatively 
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low temperature, while the confinement-deconfinement transition takes place at higher tem- 
perature as a first order transition with the discontinuous change in the value of ^. In the 
PNJL model, these two distinct transitions interfere each other through the quark quasi- 
particle loops. As a result, the two transitions takes place at a similar temperature: the 
transition temperature for deconfinement becomes lower and becomes smooth cross-over 
transition, while the chiral transition takes place at higher temperature, since the quark 
excitations are suppressed at low temperatures by small expectation value of the Polyakov 
loop. This important observation was first made by Fukushima.'^ With a finite bare quark 
mass mo, the chiral symmetry is broken explicitly and the chiral transition also become a 
smooth cross-over transition. 




Fig. 1. Temperature dependence of effective quark mass Mq and expectation value of the Polyakov 
loop: the left panel is for vanishing bare quark mass (the chiral limit) and the right panel is 
with finite bare quark mass mo which is chosen to reproduce the pion mass rriTr = 140MeV. 



15. Mesonic correlations: collective v.s. non-collective 



Now we come to our main task to evaluate mesonic correlation energy. The contribution 
of the mesonic correlation energy to free energy, ^m, or the corresponding pressure Pm = 
—Qm/V, is given by 

(5-1) 



T 6^1 
i7Af(T, A4) = — Trjv/ln ytT, 



and 



T 6^1 
PM{T,Ai) = -— TrMln -^tT, 



(5-2) 



<t>=4>o 



where the trace is over the arguments of the shifted auxiliary fields (pi with a periodic 
boundary condition along the imaginary time axis. Using (12 -jp and (I2-8P we find that the 
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mesonic pressure is given by 

+31n ((2G)-i - il^K, g, A^))] } (5-3) 

where = 2n7rT is the bosonic Matsubara frequency and 

il^K, g, A4) = 77^(^4) + i72(a;„, g, A4) (5-4) 
77,(0;,, g, A,) = nl{A^) + i7^(a;„, g, A4) (5-5) 

are the meson self-energies with 



and 



i72(w,, g, A4) = [ul + q2 + 4M2) F^, q, A4) (5-7) 
i72(a;„, g, A4) = K + q') ^K' ^4) (5-8) 



where we have defined a dimensionless common multiphcative factor by 

F(c.., q, A,) = 2TY: tre | (2vr)3 + ,^4)^ + i^^ ■ [{I + ,^4 + c.)^ + 



(5-9) 



We now show that the contributions of the mesonic correlation to the free energy fl5-ip 
or the pressure (15-21) indeed contain those of free meson gas composed of massless pions and 
massive sigma mesons. For this purpose we first eliminate the "tadpole" terms and 
by using the stationary condition, or the "gap equation" (13- 7p which determines the quark 
quasiparticle mass Mq in the symmetry broken phase. We then find, 

{2G)-' - i7,(a;„, g. A,) = {ool + q^ + 4M^)F{uJn, g, A,) + (5-10) 
(2G)-i - i7,(a;„, g. A,) = {ul + q')F{u,,, g. A,) + (5-11) 

5.1. Chiral limit: rriQ = 

We first consider the chiral limit, mo = 0. In this limiting case, 

(2G)-^ - n^itun, g, A,) = {col + q' + 4M2)i"K, g, A4) (5-12) 
(2G)-i - i7,(u;„, g. A,) = {u^ + q^)F{u^, g, A4) (5-13) 



12 



so that we can separate the contributions of the collective bare meson modes from non- 
collective individual excitations: 

Pm{T, a,) = -^Tj2l^ [In {f3'{u;l + + AM^)) + 31n {(3'{col + q')) + 41nF(a;„, q, A,)] 

(5-14) 

The first two terms contain contributions to the pressures, respectively, of free sigma meson 
gas with mass = 2Mq and of free massless pion gas with 3 isospin degeneracy. Indeed, 
one can perform each bosonic Matsubara frequency sum by contour integration and find. 



(2vr) 



T^/BK) + 3x^/s(g) 

6UJ„ 3 



(5-15) 



where oOg = ^/q^ + 4Mq is the sigma meson energy and /_b(w) = l/{e^^ — 1) is a bosonic 
single particle distribution function. The first term is the mesonic vacuum pressure p%, 

Pm = -^J ^[3g + c^,] (5-16) 

where we have indicated an additional cutoff Ah which need to be introduced to regularize 
otherwise-divergent mesonic momentum integral. This vacuum mesonic pressure needs to 
be removed by renormalization. The momentum integral for the pressure of the massless 
pion gas can be evaluated analytically and one finds a familiar Stephan-Boltzmann pressure 
with three-fold isospin degeracy: 



d^q q 



P& = y X = 3 X (547) 

Now we evaluate the correction to the free meson gas due to underlying quark substruc- 
ture of mesons. This is contained in the additional contribution, 

ApMiT,A,) = -2r^ J ^lnFK,g,A4) (5-18) 

where Un is the bosonic Matsuibara frequencies, Un = 2TinT . The function F defined by fl5-9p 
also contains the sum over the fermionic Matsubara frequencies = (2n + \)tiT. These 
sums are computed by the contour integration on a complex z plane with multiplicative 
complex function ±l/(exp(z/3) ± 1) which has poles on the imaginary 2;-axis at e„ (a;^). We 
find 

f dz f d^a 1 

where 

J^{z,q,Ai) = F{-iz,q,A^) (5-20) 



2 
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Noting that the function T has poles only on the real 2;axis, we modify the integration path 
continuously to the one which encircles the real z axis clock-wise. Rewriting the real value 
of z by cu, we obtain 

Avm[T,A^) = -2 — ^ -In- 



■oo 
oo 



27rz 7 (27r)3e^/^-l T{uj-ie,qA4] 



r ^ 2 T{u + te,q,A,) 

where we have used the fact that the function J^{uj, q, A4) is a even function of uj and 
l/(e-'^^ - 1) = 1 + l/(e^^ - 1). 

Similarly, the sum over the fermionic Matsubara frequency in fl5-9p can be performed by 
contour integration: 

r d^p r dz 1 1 

:F[u,q, 4) - y ^ J^2^^e^P + 1 [(-{z + igA.f + Ej] ■ [-{z + igA, + + 



tr^ 



d^p 



(27r)3 2Ep2Ep+g e(^p-*^4)/3 + 1 ^^^u + _ ^^^^ w - + Ep^ 

where we have performed the contour integration by modifying the path to the paths en- 
circling the poles on the real z axis, showing only one of four terms. Again, the external 
gauge fields appear as a phase factor in the quark distribution function in J^{ujn, q, A4) as 
in the mean-field approximation. We therefore replace these phase factors by the Polyakov 
loops and then substitute them by statistical average as in (14- 2p performing sum over bosonic 
Matsubara frequencies by contour integration, we found 



where 



J^sc-A^,q) = 3 



J^{u, q) = {J'iuj, q, Ai)) = J'scA^, q) + -^pair(w, q) (5-22) 
d^p 1 



(27r)3 2Ep2Ep+q \uj + Ep — Ep+q u — Ep + Ep^g^ 

X (UiEp) - MEp^q)) (5-23) 
d^p 1 



(27r)3 2Ep2Ep^q \u + Ep + Ep^q uj — Ep — -Ep+g, 

X (1 - U{Ep) - MEp^q)) . (5-24) 



The detail of this computation is given in the appendix A. We may interpret these correlation 
energy as non-collective fluctuations of the system carrying mesonic quantum numbers with 
quenched quark distributions. 

For the computation of the correlation energy or pressure, it is convenient to decompose 
the function J^{uj ± ie, q) into real part J^iiuj, q) and imaginary part ^-"2(1^, q)'- 



J^{u ± le, q) = J^iiu, q) ± iJ^2{oo, q) = ^^1(0;, g)^ + 7-2(0;, g)2e±'^("''?), (5-25) 
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where the argument is given by 

(f){uj,q) 



tan" 



-1 ^2(w,g) 



(5-26) 



The pressure arising from the non-collective or individual excitations of the system is 
given by 



Apm{T) = {Apm{T,A^)) 



(2vr) 



°° diO 
2^ 



1 + 



ol3u} _ ]_ 



2<P{oj,q). (5-27) 



5.2. Breaking the chiral symmetry with ttlq ^ 

Foregoing analysis applies only to the chiral limit (mo = 0) and symmetry broken phase. 
For non-zero value of mo, the separation of the collective modes and the individual excitations 
is not as simple as above analysis due to the term mo/(2G'Mo) in the dispersions (15- 101) and 
(15 -111) for each kind of mesons. We expect in this case that the pions, would-be Nambu- 
Goldstone modes associated with spontaneous symmetry breaking of the chiral symmetry, 
acquire non-zero mass, m^. To calculate the contribution to the thermodynamic potential 
of such modes we need to go back to the original formula (15-31) for the pressure from mesonic 
correlation. It is convenient to write 



M^{cOn,q) = i2G)-' - {n^{cOn,q,A^)) 
M^{Un,q) = {2Gy' - {n^{Un,q,A^)) 



or using the gap equation. 



M„{un, q) = (col + q' + 4:M^){F{L0n, q, A^)) + 



2GMq 



M^{u;n,q) = {ul + q^){F{un,q,Ai)) 



mo 



(5-28) 
(5-29) 

(5-30) 
(5-31) 



2GMo 

Performing the bosonic-Matsubara frequency sum by contour integration again, we obtain 



Pm{T) 



where we have written 



d^q 1 
(27r)3 2^z 7o 

X i3 In 



e/3i^ - 1 
Mt,{uj + ie, q) 
M.„{u - ie,q) 



In 



M^joj + ie,q) 
Maiuj - ie,q) 



M^/a{uj,q) = M^/^{-iuJn,q). 



(5-32) 



(5-33) 



with {F{ujn,q, A4)) in M.„/a{—iuJn,q) replaced by (l5-22p . Expressing the integrand of (I5-32|1 
in terms of the arguments of complex JliT,/„{uj ± ie, q), we have 



Pm{T) 



d?q 
(2^ 



2^ 



1 + 



_ iJ 



[3</)^(a;,g) + </'a(w,g)] 



(5-34) 
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where 



0V<^(w, = tan"^ I " (5-35) 

with 

M^,.{uj ± ze, q) = Ml/,{uj, q) ± iMl/,{uj, q). (5-36) 

We note that unhke the case of the chiral hmit, (p^^ and </)cr defined as eq. fl5-35p contain 
contributions not only from non-collective modes but also from collective modes. Here we 
call the collective meson mode only when it appears as an isolated pole on the real u axis: 
it requires the condition that both real and imaginary parts of M.T^/a{uj^q) vanish for a 
certain combination of uj and q determining the mesonic dispersion relation. In this case 
the argument 0(0;, q) becomes a 5 function, 5{(jJ — Ug) with the mesonic dispersion relation, 
LOq ^ ^/ q'^ + with m = mT^inifj). 

On the other hand, if only the real part of A^(cj,g) vanishes while the imaginary part 
is non-vanishing, such mode decays by Landau-damping into non-collective excitations. We 
refer these modes as the non- collective individual excitations As the temperature increases 
the low energy boundary of the continuum of the particle-antiparticle excitations comes 
down and eventually absorbs the meson poles. Hence no isolated meson poles remain. In 
the chiral limit this happens at T = Tc where the chiral symmetry is restored and pions, 
Nambu-Goldstone modes, disappear. In this case, there is no meson pole contribution and 
only contribution from cuts remains in the mesonic correlation energy. With explicit chiral 
symmetry breaking by the bare quark mass, we have to go back again to the original expres- 
sion (15 •3p in order to compute the mesonic correlation energy as has been done We 
found numerically that the real part of M.{(jJ, q) does not vanish at high temperatures for all 
values of (w, g), showing no sign of persistence of mesonic modes, as shown in the section 6. 

5.3. Cut-off parameters 

Since our model contains some divergences, we have to introduce cut-off parameters 
to suppress these divergences. For calculating pressure of mesonic correlation, we need 
two cut-off parameters. First divergence appears under the mean field approximation. A 
cut-off parameter is introduced in the eq. fl3-4p . This is the vacuum pressure under the 
mean field approximation. The first term of eq. fl3-4p diverges because of the momentum 
integral. We remove higher momentum modes by the cut-off Af. The thermal excitation 
term, the third term on the right hand side of eq.f l4-lip . doesn't need any cut-off because 
quark distribution functions suppress the excitations of high momentum modes. Hence, 
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when once Aj is introduced in the vacuum term, we can calculate the pressure under the 
mean field approximation without any more divergence. 

However, when we calculate pressure of the mesonic correlation, additional divergence 
appears in the momentum integral with the first term in the bracket of eq. fl5-27p . This 
divergence arises from the meson momentum integral and we cannot eliminate it by quark 
momentum cut-off Af alone. This problem arises because NJL type models are unrenormal- 
ize: even if we suppress a divergence from a calculation of some order, additional divergence 
appears in the next order calculation. For this reason, we have to introduce second cut- 
off parameter Ab to suppress this divergence from bosonic loop calculation. Since Ab has 
no relation to Af mathematically, we should determine the value of it based on physical 
considerations. In this work, we choose Ab = 2Af^ 

The integral with the second term in the bracket of eq. fl5-27p doesn't diverge. Since 
bosonic distribution function as a function of u suppresses excitation modes at large u 
region, no divergence appears in the u integral. As for q integral, it seems at first sight 
to diverge because the bosonic distribution function is only a function of u. However since 
becomes very small at large q, the divergence doesn't appear. One can see this from 
eq. f lB-lll) . When we fix u and see large q region, it is always space-like region so that 
the behavior of J-2 corresponds with J^scatt,2- From eq. (IB-lip . we see that Nscati.^) becomes 
exponentially small at large q region. As a result J^scat,2 also gets close to zero from eq. flB-16p . 
Although eq.f lB^Tj) is an approximation under the condition of small e, this approximation 



is justified as long as q becomes large for fixed u in the q integral of eq.f l5-27l) . 

§6. Numerical results 

6.1. Pressure 

We show in figure [2] the pressure calculated by the present method. The left panel is 
the result computed in the chiral limit mo = 0. At low temperatures below the second 
order chiral transition temperature Tc, the pressure from quark quasiparticle excitations are 
suppressed by the quenching of the distribution function: only excitations of triad of quarks 
with color singlet configuration are allowed with effective excitation energy three times that 
of massive quark excitation, this contribution to the pressure is hence suppressed. The 
dominant mode of excitations which determine the pressure is due to massless pion excitation 
with 3 isospin degeneracy as shown in the red solid curve. We found the contribution from 
massive sigma meson excitations and non-collective individual excitations of the quark triads 
and antiquark triads are also negligible at low temperature. At high temperatures, massless 
quark excitations in the mean-field approximation dominate the pressure with comparable 
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mean-field (quark+gluon) 






mean-field (quark only) 





















meantield + meson correlation ■ 
mean-field (quark-i-gluon) ■ 
mean-field (quark only) ■ 




Fig. 2. Pressure scaled by as a function of temperature: the left panel is for vanishing bare 
quark mass (the chiral limit) and the right panel is with finite bare quark mass rriQ which is 
chosen to reproduce the pion mass = 140MeV. In comparsison, the pressures calculated in 
the mean-field approximation are shown with (without) gluon contributions. 



contributions of the gluon pressure phenomenologically introduced in the construction of the 
effective potential. Contribution from the mesonic correlations becomes negligible at high 
temperatures, and the pressure approaches to that of ideal gas of massless quarks and gluons 
in our model. We note that at large temperatures the cut-off A in the quark momentum 
integration reduces the pressure compared to that of free quark-gluon gas, as seen in the 
effective reduction of the Stephan-Boltzmann constant p/T^ from a = 3x2x2x2x7/8x 
7r^/90 = 2.3. With symmetry breaking finite bare quark mass, the pion becomes massive 
and its temperature decrease exponentially at very low temperature, otherwise the behavior 
of pressure at higher temperatures is qualitatively unchanged. 

6.2. Collective modes and non- collective modes 

Since mesonic excitations contain two types of modes, collective meson modes and non- 
collective individual excitations, we consider in this subsection the two contributions sepa- 
rately, for the chiral limit where the separation is manifest, and for the case with the chiral 
symmetry breaking, where we need to separate them carefully due to the mixing via the 
bare quark mass term. 

6.2.1. Chiral limit 

As shown in figure 2, the pressure at low temperatures satisfies the Stephan-Boltzmann 
relation p/T^ ~ 0.33. It means that pressure at low temperature comes only from pion 
collective modes, eq. fl5-17l) . and the contributions from sigma meson collective modes and 
non-collective modes are very small. The contribution from thermal excitations of sigma 
meson is suppressed compared to that of pions due to large sigma meson mass, m„ = 2Mq. 
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The contribution from non-collective modes is also expected to be suppressed due to the 
same reason: the effective mass of quark triad is 3Mq as it appears in the quenched quark 
distribution. Since the contribution of non-collective modes to the pressure comes from the 
phase space integral of the argument (f){u, q) of J^{u + ie, q), we like to inspect in detail the 
behavior of J^{uj + ie, q). 




Fig. 3. Kinematical regions where J- 2 is non- vanishing is shown by shaded areas: (1) < q 
(scattering), (2) w > \/ + (2Mq ) (pair excitation). The left panel is for T = 0.7Tc, while the 
right panel is for T = lATc- 

We show first the region where the imaginary part of function J^{u + ie,q), J^2, be- 
comes non- zero in figure [3l It consists of two regions, one in the space-like region (w < q) 
corresponding to the scattering, and the other in the time-like region {u > q) where pair 
excitations contributes. The upper boundary of scattering region is at the light cone [uj = q) 
and the lower boundary of pair excitation region is given by lo = \/ q'^ + (2Mo)^ where q 
is the momentum of the mesonic excitation and Mq is the constituent quark mass deter- 
mined by the gap equation fl3-8p . In the region between the two boundaries, J^2 becomes 
zero. There are no physically allowable excitations in the kinematical region between these 
boundaries, besides meson poles which are not shown here. As temperature increases, the 
upper boundary goes down because increase of temperature causes decrease of constituent 
quark mass Mq as a result of restoration of the chiral symmetry. Hence there appears indi- 
vidual excitations filling this region at high temperature as shown in the right panel of figure 

El 

Next we display J^i and J^2 as a function of u/q for a fixed value of g = /!/ at two 
characteristic temperatures in figure H] (T = 0.7Tc) and figure [5] (T = lATc). Below Tc, 
crosses the horizontal axis and there are two valleys aX uj/q = 1 and w/g ^ 1.6, on the other 
hand, above it doesn't cross the horizontal axis and the valley at w/g ^ 1.6 goes down 
and is absorbed into continuum. 
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Fig. 4. Left: T\ as a function of w/g at T = 0.7Tc. Right: J-2 as a function of w/g at T = 0.7Tc 




Fig. 5. Left: T\ as a function of w/g at T = 1.4Tc. Right: J-2 as a function of w/g at T = lATc. 

In order to see the behavior of (f){u}, q), we plot Argand diagrams of J^{uj±ie, q) in figure O 
Integral of (f){oj, q) determines contribution to pressure of non collective modes. As we have 
seen in subsection I5.3[ q integral of </> doesn't have any divergence at large q. One can see 
in figure [6] that J-2 becomes small as q/u becomes large with fixed u at both temperatures; 
actually J^2 and q) become exponentially small at large q. 

In figure Elthe point marked by (1) is in the time-like region {u > q) and points marked 
by (3) or by larger numbers in the parentheses are in the space-like region {u < q). /^From 
upper panel, we can see that (p at low temperatures becomes almost vr in the space-like region 
because is given by eq. fl5-26p . It seems that the contribution to pressure from the time- 
like region is large. However, J^2 in the time-like region contains large vacuum contribution 
much larger than that of thermal excitations. The vacuum pressure will be subtracted so 
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-0,01 I 1 1 1 1 1 

(I 0115 III (1,1 5 11,: (1.25 

Fl 

Fig. 6. Argand plots of {Ti{ijj,q).J'2{^,Q))- The upper panel is calculated at T = 0.7Tc, while the 
lower panel is at T = lAT^, each with 3 fixed values of co indicated by different colors. Numbers 
in the parenthesis correspond to 6 selected values of q/oj at (1) 0.8, (2) 1, (3) 2, (4) 5, (5) 10 
and (6) 20. 

that the contribution from the time-like region becomes small. Namely, the effect from finite 
temperature appears in the space-like region. Moreover, the value of J^2 is much smaller 
than that of Ti. As a result, the value of itself becomes very small. It means that the 
contribution from non- collective modes is negligible. At high temperature, the value of J^2 
is lager than that at low tempe ratures because constituent quark mass becomes zero due 
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to the restoration of the chiral symmetry and because quarks which are suppressed by the 
Polyakov loop at low temperatures become to excite at high temperatures. However the 
value is smaller than that of T\ so that the value of is small. In addition, a cancelation 
between positive value of and negative value of occurs. As a result, the contribution 
from non-collective modes at hight temperatures is also small. 

6.2.2. The case with mo 7^ 

In the case of the chiral limit, we find that the contribution from the collective modes 
among mesonic excitations are much larger than that of the non-collective modes. We antic- 
ipate that it remains unchanged in the case which quarks have finite bare mass, violating the 
chiral symmetry explicitly. In this subsection, we show that even if quarks have finite mass, 
collective modes also exist at low temperatures and they disappear at high temperatures. In 
addition, we show the largest contribution at low temperatures still comes from the collective 
modes and the contribution from the non- collective modes remains negligible by calculating 
integral of argument of complex function A^(ciJ ± ze, q). 

When quarks have a finite mass, we cannot separate collective modes from non-collective 
modes easily by factorization as in eqs. fl5-12l) and fl5-13p . In this case, we first need to 
determine whether collective modes exist or not. For this reason we calculate the real and 
imaginary part of M.{uj ^ i^,<l) separately given by eq. (15-361) . The conditions of isolated 
meson poles are given by the vanishing both A^i and A^2 ; in particular the condition 
Aii{uj,q) = determines the dispersion relation of collective meson modes. 

We show the plots of Aii and A^2 of pion as a function of u scaled by q in figure [71 At 
low temperatures, Aii crosses the horizontal axis in the time-like region, and A^2 is also zero 
at the same value of u/q where A4i vanishes. It means that a collective pion mode exists 
at this point. Furthermore, there are no individual excitations in the region where = 0, 
since A^2 is proportional to ^"2(0;, q), so that the collective pion modes appear in this region 
as stable excitations with infinite lifetime. On the other hand, at high temperatures, the 
region where A^2 = disappears. It means there are no isolated pion poles with infinite 
lifetime at high temperatures even though Aii becomes zero. We also show Aii and A^2 of 
sigma meson as a function of u scaled by q in figure [SI The situation is same as the case 
of pion. At low temperatures there is an isolated sigma meson pole in the time-like region, 
while there are no poles at high temperatures. 

Next, we consider the contribution from the non- collective modes. In the same way as the 
case in the chiral limit, it is useful to see Argand diagrams of Ai{uj±ie, q) for understanding 
the contribution to the pressure from non-collective modes. Note that unlike the case of 
chiral limit A^(a; ± ie,q) contains not only non-collective modes but also collective modes. 
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However, since collective modes always exist in the region where A^2 is zero, we can separate 
contribution of the collective meson excitations from that of the non-collective modes. The 
pressure from collective meson excitations is given by that of ideal gas with the modified 
dispersion relation determined by Aii{u, q) = 0. Therefore the pressure calculated from the 
integral of (p defined by eq.f l5-35p corresponds to the pressure of the non-collective modes. 

Argand diagrams of AiT,{u + ie, q) and A^o-(a; + ze, q) are shown in figure [9] and figure 
[TD] respectively. From these figures, one can see that M.2 of pion and sigma meson becomes 
small as q increases. Actually becomes exponentially small at large q for a fixed u so 
that the q integral converges. The behavior of M..„{(jJ + ie,q) and A^o-(i^ + i^,q) at high 
temperatures is the same since Mq = when the chiral symmetry is restored. 
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Fig. 8. Ml and M2 of sigma meson as a function of lo scaled by q.The upper panels are calculated 
at r = 0.7Tc and the lower panels are calculated at T = lAT^. 

§7. Conclusion 

We have studied the quark-hadron phase transition at zero net baryon density in terms 
of the PNJL model which is designed to incorporate both the chiral phase transition and the 
deconfining transition expected from QCD at finite temperature. We have shown that the 
model describes the low temperature phase as a gas of light mesons and the high temperature 
phase as a gas of free quarks and anti-quarks with restored chiral symmetry. The transition 
between two phases is continuous cross-over transition and the thermally excited degrees of 
freedom are transformed from mesonic degrees of freedom to the quark degrees of freedom 
continuously as temperature increases. 

The quark degrees of freedom is suppressed at low temperatures due to the phase cancel- 
lation of the color gauge fields acting on different color modes of quark excitations, turning 
the quark distribution functions to the distribution of the quark triads whose contribution 
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Fig. 9. Argand diagrams of A4Tr{uj + i€,q). The upper panel is computed at T = 0.7Tc and the 
lower panel is cmputed at T = lATc- 



to the pressure is suppressed exponentially due to the large effective mass. The pressure of 
the low temperature phase is dominated by the mesonic excitations which arise as collective 
mode of excitations of the mean thermal distribution of the quark triads; the contribution of 
the non-collective excitations of quark triads are suppressed. As the temperature increases 
the collective meson modes are absorbed into the continuum of individual excitations. They 
may persist as resonances with finite lifetime due to the Landau-damping. We found the 
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Fig. 10. Argand diagrams of + ie,q).Tlie upper panel is computed at T = 0.7Tc and the 

lower panel is cmputed at T = lAT^. 

contribution of these modes to the pressure is neghgible at high temperatures as welL In 
the high temperature deconfining phase, our model gives essentially the same results as the 
Nambu-Jona-Lasinio model: the equation of state becomes that of non-interacting massless 
quark-gluon gas. 

We note that the quark triad is still not a baryon. There is no spatial correlation of 
three quarks since they share exactly the same amount of momentum and energy. They also 
have somewhat over-reduced degrees of freedom. It is important to remedy these points for 
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studying further implications of the model at finite baryon chemical potentials. 



Acknowledgements 

We thank G. Baym, J. P. Blaizot, H. Fujii, K. Fukushima, T. Hatsuda, F. Lenz and B. 
Svetitsky for useful discussions and comments on this work. In particular, we thank G. 
Baym and T. Hatsuda for illuminating discussions on the effective degrees of freedom of 
quark triads. KY's work has been supported by the Global COE Program "the Physical 
Sciences Frontier", MEXT, Japan. TM's work has been supported by the Grant-in- Aid # 
21540257 of MEXT, Japan. 



Appendix A 

Pressure by non- collective mesonic excitations in the chiral limit 



Here we present the detail of the computation of the pressure exerted by the non-collective 
mesonic excitation modes given by (5.18). 

We first perform the sum over the bosonic Matsubara frequencies a;„ in (15-181) by contour 
integration in the complex z-plane: 

d^q 1 r , 1 



Apm{T,A^) = -2T 



dz- 



-\n.F{—iz^ q, A4) 



(Al) 



(27r)3 27ri Jc e^^ - 1" 

where the contour C encircles the imaginary axis counter-clock-wise (Fig. A-1). The function 
1 / (e^^ — 1 ) has poles at z = 2iinTi = Uni and is analytic everywhere else. We change the 
integration path along the contour C to the path which encircles real z-axis clock-wise (See 
Fig. A-2). Denoting the integration variable z by a; on the real z-axis and then converting 
the integration along the paths sandwiching negative u region to the integration around 
positive real axis, we find 



ApMiT.A^) = -2 
where we have written 



d^q 1 



doj 



- 1 



In 



J^{(jj + ie, q, A4) 
— ie, q, Ai) 



(A-2) 



J^{u, q, A4) = F{-iu, q, A4) 



tr. 



d^p 



1 



(A-3) 



(27r)3 [(e„, + ^^4)2 + E^] - [(e„/ + gA^ - zoo) 

where e„' = {2n' + l)7rT is the fermionic Matsubara frequencies for quark-quasiparticles. If 
necessary, we may rewrite the above expression by the partial integration, to 

'^{10 + ie, q, A4) 



ApMiT^A^) = -2T 



— / rfw w + -ln(l - e^"") — In 

{2Txf2TiijQ L p y L w - ie, g, J 



(A-4) 
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We now express the sum over the fermionic Matsubara frequencies e^' by contour inte- 
gration: 

-7-/ A \ I ^^P f dz -1 1 

J^yu, q, A4) = tTc 



(27r)3 2m e^^ + 1 [{-{z + igA^f + E^] . [_(^ + ,gA, + + E^^J 

(A-5) 

where the function ^^i^ is chosen so that it gives poles at the fermionic Matsubara frequen- 
cies on the imaginary axis; minus sign is necessary for insuring the right sign of the residues 
at these poles. The integration path is taken as the same as the bosonic contour integration. 
We change the integration path again to the path enclosing the real z-axis and perform the 
integral around the poles of the integrand at the real values of z. The result may be written 
as a sum of two terms: 

J^{U, q, A4) = j;cat(w, q, A4) + Jpair(w, q, A4) (A-6) 

the first term is scattering term 



uj + Ep — Ep+q CO — Ep + Ep^ 
xtr, {f{Ep - igA^) - f{Ep+, - igA^)) (A-7) 

while the second term corresponds to the pair creation and annihilation term 

d^p I { I 1 



. , f d^p 1 / 



CO + Ep + Ep^g U — Ep — Ep+g 



xtr, (1 - f{Ep - igA^) - f{Ep+g - igA^)) (A-8) 

Again, the external gauge fields appear as a phase factor in the quark distribution function 
as in the mean-field approximation. We therefore replace these phase factors by the Polyakov 
loops and then substitute them by statistical average as in (14- 2p : 



J^{u;, q) = q, A4)) = J^^cA^, q) + ^pair(c^, q) (A-9) 



where 



^pair(w, g) = 3 



d^p 1 f 1 1 



(27r)3 2Ep2Ep+g \uj + Ep — Ep+g u — Ep + -Ep+g, 

X (UiEp) - MEp^g)) (A-IO) 
d^p 1 f 1 1 



(27r)3 2Ep2Ep^g \u + Ep + Ep^g uj — Ep — -Ep+g, 

x{l-U{Ep)-U{Ep+g)) (A-11) 
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We may interpret these correlation energy as non-collective fluctuations of the system car- 
rying mesonic quantum numbers with quenched quark distributions. 

For the computation of the correlation energy or pressure, it is convenient to decompose 
the function J^iuj ± ie, q) into real part J^i{uj, q) and imaginary part J^2{^, q)'- 



J^iu ± te, q) = q) ± ij^^i^^, q) = ^ J^i{u:,qY + J^2{u:,qYe^"^^^''^\ (A- 12) 

where the argument (f) is given by 

0(a;,g) = tan-i^^^ . (A- 13) 

The real part and imaginary part of the function J-" are further decomposed into two parts: 
the scattering term and the pair excitation term. The two components of the real part are 
given by the principal part integrals: 

f d^p 1/1 1 \ 

J {2iry 2Ep2Ep+g \u + Ep - Ep+g u-Ep + Ep+gJ 



X iUiEp) - U{Ep+g)) 
d^p 1 f 1 1 



(27r) 3 2Ep2Ep+g \u + Ep + Ep+g uj-Ep- Ep+g , 

x{l-U{Ep)-U{Ep+g)) (A-15) 

while the two components of the imaginary part contains the energy conserving 5-functions: 

d^p 1 



f d-^p 1 
:F.....,M = - y i2.f2Ep2Ep^g - 



^{5{u + Ep- Ep+g) -6{uj-Ep + Ep+g)) 
f d^p 1 

X {6{u + Ep + Ep+g) - 5{u -Ep- Ep+g)) (A- 16) 

Similar function for the finite temperature bosonic system has been computed.'^S'ESJ It 
is evident that the scattering term has non-zero imaginary term in the space-like energy- 
momentum region (q^ > cj^), while the pair creation/annihilation term is non- vanishing 
only in the time-like region (q^ < u"^). It is important to note that non-collective mesonic 
correlation arises only from non- vanishing imaginary part of J^{q,oj). 

Appendix B 

Computation of J^sca.t{^ , l) o-nd J^p^ni^jJ , q) 



Here we present the computation of the phase space integral in J^sca.t{^, l) and J-pair(i^, q) 
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We first adopt the polar coordinate system, d^p = dpp^ d cos 6 d(f), and then change the 
variables (p, cos^) to new variables {E,e) defined by 

E = ^{E, + E,+^) (B-1) 

^ = Ep+q — Ep. (B-2) 

Note that 2E represents the energy absorbed by pair excitation while e is the energy trans- 
fered by scattering. The Jacobian of this transformation is 

d{E,e) ^ p\ 
d{p, cos 6) EpEp+q 

so that the integration measure absorbs the Lorentz factor EpEpj^q. 

d^'P dpp'^d cos Odd) 1 ,, 

—f- = ^ = -dEdedct> (3-4) 

Note that by construction 

2Ee = El^^ -Ep = 2pq cos + (B-5) 
so that the domain of the integration on the {E, e) plane is constrained by the condition 

Ee-^ 

cos^d = ^<1 (B-6) 

pq 



Since p = ^E"^ — = y (i? — — m'^, this implies 

E^e^ < g2 (^E^ + _ ^2 _ ^3.7) 

The region of the integration on {E, e) plane is indicated by gray zone in Fig. [TTJ 

We can now rewrite (lA-lip and flA-ll|) in simpler, physically more transparent forms. 
The scattering term flA-lip may be written as 

^scat(w, g) = 7^ [' de (— —] iVseat(£) (B'S) 

where 

iV,eat(e, q) = ^j^dE (^U{E - ^e) - U{E + ^e)^ (B-9) 

may be interpreted as the number of scattering states with energy- momentum transfer {e, q). 
The lower boundary of the ii^-integration, Eq, is given from the constraint flB-7|) as a function 
of (e,g): 



i?o(^,g) = fJl + ^ (B-10) 
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Fig. 11. Region of integration on (e, E) plane is shown by gray area. 



If the condition e << q is fulfilled, the integration may be performed approximately as 

(B-11^ 



A^scat(£,g) ^ -7; 

27r q 



dE 2tt q 



Note that the number of scattering final states decreases exponentially for large q. Similarly 
we can express the pair excitation term (lA-lip as 

1 /■ .„ / 1 1 



-7^pair(w, q) 
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dE 



Eq 



u~2E U + 2E 



Np^iriE, q) 



(B-12) 



where 



Np.Sr{E, g) = ^ de (^1 - ME - ^s) - U{E + ^e)^ (B-13) 

may be interpreted as the number of pair excitation states with energy- momentum (2i?, q). 
The boundary of the ^-integration is obtained from flB-7p as a function of E and q: 



e±{E,q) = ±qJl 



(B-14) 



E^ - q^/A 

Since the value of e is constrained by —q < e < q, we can estimate the integration approxi- 
mately and find 



iVpair(^,g) 



1 'e^-e^){l + 2U{E)) = -Jl 



2'Kq 



E^-q^A 



[l-2Um (B.15) 



The functions A''scat(£^) and Npair{E) determine the imaginary part of J^sca.t{^ + q) and 
•^pair(i^ + q), respectively. Using the expression (IB-Sp . we find 

1 n 

-7%,cat,2(c^, 1) = ^ J (-vr5(a; - e) + 7iS{u + e)) A^scatl^) = -A^scat(w) (B-16) 
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where we have used that A'^gcatlf^) defined by (]B-9p is an odd function of e when extended to 
the negative value of e. Similarly, from (1B-12P we obtain 

•^pair,2(c^, q) = ^ [ dE {-7r6iu - 2E) + Tx5{uj + 2E)) N^^UE) = -iVpair(c^/2) (B-17) 

Appendix C 

Separation of collective mesonic modes with finite bare quark mass 



In the presence of non-vanishing bare quark mass, the collective meson modes are also 
contained in fl5-32p or 



Apm{T, a,) 



(Pq 1 
{2Tif2^i Jo 



duj 



1 + 



X <^ 31n 



A^^(w + ie, q) 
M-„{u - ie,q) 



+ In 



Majco + ie,q) 
Maiuj -ie,q) 



(C-1) 



together with non-collective excitation modes. In this case we can still separate the contri- 
bution from the collective meson modes. 
We first note that for small, but finite e. 



Ml{u, q) = i-u^ + q') q) + 



_mo_ 
2GMo 



Mliu, q) = i-u' + q' + AMl)J^,{u, q) 



2GMq 



while 



(C-2) 
(C-3) 



Ml{u, q) = -2ujeTx{uj, q) + (-w^ + q^) J'2(w, q) 



Ml{oJ,q) = -2ujeJ^i{uj,q) + 



-UJ 



q2 + 4M2)J-2(a;,g) 



(C-4) 
(C-5) 



so that in the kinematical region where no individual excitation exists, J^2{^, O) = 0, we have 



Ml,AuJ,q) = -2ujeJ^i{uj,q) 



(C-6) 



Since A^^/^ = at meson poles, the argument of A^o-(i^ + i^,q) becomes ±7r/2, depending 
only on the sign of J^i(cj, q). In the limit e — )■ 0, the poles move to the real u axis. 
The locations of the meson poles on (cj, q) plane are determined by the conditions: 



Mlioo, q) = {-uj^ + q2 + 4Mo') q) 



2GMp 



Mliu,q) = {-u^ + q^)Ti{u,q) 



mp 
2GMp 



(C-7) 
(C-8) 
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for the a meson pole and the pion poles respectively. The bare quark mass parameter mo 
is determined so that it generates the physical pions mass in the vacuum, hence it should 
satisfy near g = the on-the-mass-shell condition uj = , hence 

mo 



2GMn 







(C-9) 



where J-'^ad^, q) is the vacuum piece of the function F{uj, q): 
Subtracting ( ]C-9|) from (1C-7P and (IC-Sp we obtain 



M^{co, q) = i-u^ + q2 + 4M^)J^{u, q) + m^ j;,,(m,, 0) 
M^{u, q) = (-cj^ + q^)J'{uj, q) + mlJ'^ac{m^, 0) 



(C-10) 



(C-11) 
(C-12) 



The contribution of these "pole terms" to the cj- integral is essentially the same as in the 
chiral limit, except that the masses of mesons are shifted. 

The contribution from the cuts comes from the non-vanishing imaginary part of J-'{uj ± 
ie, q), which generates non- vanishing imaginary parts in A^o-(ij-' ± ^e, q) by 



Mliu, q) = {-u^ + q2 + 4M2) J-2(a;, q) 
Mliu,q) = i-u' + q^)T2{oo,q) 



(C-13) 
(C-14) 



while the real parts of JliT,/„{uj ± ie, g) remain non-vanishing on the cuts, so that we have 
the phases 



<^7r(i^, q) = tan" 



q) = tan" 



-1 Ml{u,q) 
MUu^q) 



^i Ml{u;,q) 
MUoo,q) 



tan 



J-2 



^1 



mo 



tan 



2GMo q2 - ^2 y 

J'2 



\ 



^1 



mo 



(C-15) 



(C-16) 



2GMo q2 + 4M2 - 



shifted slightly from the values in the chiral limit. These phases may be interpreted as the 
phase-shift of the scatterings of quark against quark or anti-quark.'^''23 
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